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We present a microscopic model of the quantum paraelectric-ferroelectric phase transition with a 
focus on the influence of coupled fluctuating phonon modes. These may drive the continuous phase 
transition first order through a metaelectric transition and furthermore stimulate the emergence of 
a textured phase that preempts the transition. We discuss two further consequences of fluctuations, 
firstly for the heat capacity, and secondly we show that the inverse paraelectric susceptibility displays 
~ T^ quantum critical behavior, and can also adopt a characteristic minimum with temperature. 
Finally, we discuss the observable consequences of our results. 
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I. INTRODUCTION 

Ferroelectric materials feature in many modern day 
electronic devices including computer memory and ca- 
pacitors, and are a simple setting for studying quantum 
criticality ^^^^ paper we focus on the family of 

displacive ferroelectrics where the optical lattice modes 
condense, forming a structural distortion. Near to quan- 
tum criticality excitations can become highly degenerate 
and new phases can emerge. Motivated by recent ex- 
periments that signal the emergence of novel quantum 
critical behavior in ferroelectrics we explore the pos- 
sibility that transverse components of polar fluctuating 
phonons conspire to drive a first order displacive meta- 
electric transition and investigate the implications for the 
inverse susceptibility. 

The soft-mode optical phonons in ferroelectrics can be 
well-described by a bosonic field theory. If the dynamics 
were not damped by free electrons and the interactions 
remain short-ranged then the general quantum critical 
behavior would adhere to the well-established rules re- 
viewed in Ref. However, in ferroelectrics the motion 
of the atoms in optical modes leads to the emergence 
of electric dipoles. A good description of these long- 
range dipole forces is essential to properly describe the 
ferroelectric transition. The effect of long-range dipolar 
forces was first studied by Rechester Q and Khmel'nitskii 
and Shneerson Aharony and Fisher Q found that 
anisotropics associated with the dipolar interaction led 
to a universality class in the classical ferroelectric. The 
quantum ferroelectric phase transition in the mean-field 
approximation, and its universality class, was studied 
by Roussev and Millis Q. However, recent experimen- 
tal evidence points to new physics that emerges close 
to quantum criticality; for example the coexistence of a 
quantum paraelectric phase with a quantum ferroelectric 
phase in ^^0-exchanged SrTiOs provides strong evidence 
for a first order phase transition [l| . Additional motiva- 



tion to study ferroelectrics arises from the inverse dielec- 
tric constant behavior of SrTiOa which falls at low tem- 
perature before increasing as ~ at intermediate 
temperatures and rises as ~ T at high temperature. 
One suggestion is that new phenomena are driven by the 
couphng of acoustic to optical phonons However, 
inspired by the ramifications of quantum fluctuations in 
ferromagnets [lol |. we show that the transverse coupling 
of fluctuating phonons can drive a first order metaelectric 
transition. 

Having realized that fluctuations can cause the emer- 
gence of a first order transition it is natural to search for 
further phase reconstruction. Motivated by the develop- 
ment of a textured FFLO phase [ll|, , and evidence for 
a textured ferromagnetic state near to the ferromagnetic 
first order transition [3, Ell i here we search for the 
emergence of an analogous textured ferroelectric phase. 
Finally, to connect to prevailing experimental methods, 
we derive an appropriate expression for the inverse sus- 
ceptibility that is consistent with recent experimental re- 
sults 0, [H, over a wide range of the phase diagram, 
and demonstrate that the transverse coupling of fluctuat- 
ing phonons could cause it to have a characteristic mini- 
mum at low temperature. 
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II. ACTION AND MEAN-FIELD THEORY 



We adopt a bosonic field theory to describe the soft 
optical phonon modes that should recover the main phys- 
ical behavior of the system. The order parameter of the 
theory is the local polarization <^(x, t) = "^2^=1 GiVi{'x., t), 
which is formally defined for one unit cell at x containing 
n atoms of charge each individually displaced through 
Vi by the optic mode. As the optical phonon softens, the 
action develops an instability and the order parameter 
must describe both thermal and quantum fluctuations. 
Following Roussev and Millis Q we describe the action 
in three-dimensional space and imaginary time via the 
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FIG. 1: (Color online) The phase diagram in the u ~ v plane 
at zero temperature in the mean-field approximation. The 
cross-hatched forbidden region denotes where the polarizabil- 
ity would diverge without higher order corrections. The solid 
thick line highlights a first order phase boundary between the 
light grey region that denotes diagonal order, and the dark 
grey which labels the Ising phase. In each regime the inset 
axes illustrate the polarization solution highlighted by the red 
vector. 



Ginzburg-Landau phenomenology 
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where a is the lattice constant, c is the speed of the 
phonons, = q^, the dimensionless momenta — tt < 
qa < TT, and the second summation is carried out under 
the conservation of momentum (qi + q2 + qa -I- q4 = 0). 
Since the field <p describes an electric dipole, the action 
includes a long-range dipole interaction, and also a cou- 
pling to the underlying lattice through the parameters r, 
/, g, and h. The terms u and v that describe the local 
anharmonic interactions give a net positive contribution 
which ensures that the polarization remains bounded. In 
general these parameters are tensorial, but for simplic- 
ity we have assumed that they adopt cubic symmetry. 
Estimates for the parameters shown in Table |T] were ob- 
tained from ab initio calculations [1, [TM Tot in the two 
key ferroelectrics SrTiOa and KTaOa Q]. The typical 
energy scale of ferroelectric fluctuations along (100) is 
Eq = hnc/a] using this definition we can then employ a 
dimensionless bosonic Matsubara frequency Cj — oj/Eq, 
and a dimensionless temperature T = T/Eq. Through- 
out the paper we adopt the units a = h = = 1. 

To establish the connection to previous work we first 
consider the mean-field phase diagram that is sketched 
in Fig. [TJ Making the ansatz that the ground state is 
uniform we obtain the action S = rcfy^ + {u + v)(p'^, where 
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arameters for the ferroelectrics SrTiOa and 



(/> = |0|. When V < the polarization (f>x — cfiy — 0, 0^ = 
—r/2(u + v) has an Ising configuration, whereas when 
V > the polarization 0^ = 0^ = (/)^ = —r/2{3u + v) 
exhibits diagonal order. The term proportional to v con- 
trols the polarization direction in the ferroelectric phase, 
whereas the u term is rotationally invariant. We note 
that whilst sweeping v through v ^ with u > the 
first order rotation of polarization direction is accompa- 
nied with a continuous change in the magnitude of the 
polarization. This is driven by a similar mechanism to 
the Blume-Emcry-Griffiths model involving two bosonic 
fields [l^l . Within the mean-field approximation the con- 
dition for stability of the polarization is that the net co- 
efficient of the quartic term is positive which translates 
to u + v > when V < and u + v/3 > Q if v > 0. If these 
conditions arc not fulfilled then higher order terms must 
be included and rather than undergo a second order tran- 
sition at r = 0, the system might have a first order ferro- 
electric transition at mean-field level. We can neglect the 
higher order terms such as At/)^ provided that the model 
remains stable, which requires that A(^^ u + v. Here 
we wish to investigate whether near criticality the fluc- 
tuating modes can conspire to drive an otherwise second 
order transition to become first order. In order to access 
this behavior we now go beyond mean-field and consider 
the consequences of quantum fluctuations on the system. 



III. FIELD INTEGRAL FORMULATION 

To account for fluctuation corrections to the system 
Roussev and Millis ^ employed the renormalization 
group, which is tailored to study the well-established sec- 
ond order ferroelectric transition. However, motivated 
by recent experiments [l|, Q we wish to explore the pos- 
sibility of a first order metaelectric transition. There- 
fore, rather than considering just the corrections due to 
slow fluctuations that arc encompassed by renormaliza- 
tion group, we need to consider fluctuations ip over all 
length scales in the polarization 4>+xp around the saddle- 
point solution c(). When u we can neglect fluctua- 
tions in "0 beyond second order which reduces the action 
to 
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(a) u — V plane, r < 



(b) u — r plane, v > 
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FIG. 2: The phase diagram at T = in the (a) u — v plane with r < 0, and (b) u — r plane with v > 0, both at zero temperature. 
The cross-hatched forbidden region denotes where the polarizability would diverge without higher order corrections, the light 
grey denotes diagonal ordered polarization, and the dark grey the (a) Ising phase and (b) paraelectric phase. Solid thick lines 
denote first order phase boundaries, dashed lines second order transitions, and the circle the tricritical point. 



where = ^Sa^p + Ua^f3, the diagonal inverse Green 
function takes the form = + + r + fq^ + (ff — 
hq'^)q'^/ q"^ + {4u+6v)(j)a+'2u4>'^ , and the off-diagonal terms 
aieUa^p = {g — hq'^)qaq0/q'^ +4u(pa4>p. We now integrate 
over quantum fluctuations to yield the free energy 



F = 



ixrlnl 
2/3 



(3) 



where f3 — 1 /T is the dimensionless inverse temperature. 
If = and r ^ g ~ /itt^ or if ^ and r ^ tt^ then 
UG 1 . In this regime we can expand the inverse Green 
function in its off-diagonal terms U using TrlnG"^ = 
TrlnG"^+Trln(l-|-GJ7) which enables us to describe the 
renormalization of fluctuations by off-diagonal coupling. 
This yields 
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where ^ [q^ + r + fql + (4w -|- 6v)(j)l + 2w<^2ji/2^ 
remove the fluctuations of the static uniform component 
of xp, which are included in 0, we must introduce the 
second logarithm. This has the effect of regularizing the 
divergence which would otherwise develop from the first 
logarithm. This expression, except for the final fiuctu- 
ation correction term, agrees with that of Ref. Q, and 
is analogous to the coupling of transverse ferromagnetic 
fluctuations that led the emergence of flrst order behav- 
ior [l3| . The condition for stability is the same as for the 
mean-field case. 

The momentum integrals are in general evaluated nu- 
merically. However, to further investigate the diago- 
nal ordered phase we make the approximation that the 
cuboid Brillouin zone boundary (— tt < qa < tt) that 
bounds the momentum space integral can be replaced 
with a spherical boundary that encloses the same total 



phase space, so has radius go = Btt^. In the low tem- 
perature limit with the polarization aligned in the (1,1,1) 
direction, the resulting integrals can then be evaluated 
analytically to yield 




r + 2{u + v)(p'^ + 4:U(j)'^/3, and were found to 
be in good agreement with the corresponding numerical 
result. 



A. Phase behavior and heat capacity 

The phase behavior of the system is shown in Fig. [21 
The forbidden region indicates where the action polar- 
izability and free energy would diverge without con- 
sidering higher order corrections to the original action. 
When considered within the framework of mean-field 
phenomenology, here the system could undergo a flrst 
order paraelectric-ferroelectric transition. However, the 
corrections due to quantum fluctuations renormalize the 
action, causing a metaelectric boundary to peel away 
from the first order transition associated with the for- 
bidden region. This metaelectric transition is consistent 
with recent experimental evidence for a first order phase 
transition [ij in ^^0-exchanged SrTiOs. In both of the 
planes considered, the line of flrst order metaelectric tran- 
sitions covers an extensive region of the phase diagram, 
terminating in a tricritical point at m = 0. The flrst 
order transition at small u is destroyed at non-zero tem- 
perature, with the tricritical point moving up the line 
of transitions to u w dTrT^^^v/r. This critical behavior 
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does not depend on the long-range dipole interactions 
since the lowest order term in o and h averages to zero 
on integrating over momenta [8|. The cj) —(j> sym- 
metry could be destroyed by applying a uniaxial electric 
field misaligned to the lattice. 

A further ramification of the quantum fluctuation cor- 
rections is that the rotation of the polarization from Ising 
to diagonal order no longer occurs where v turns nega- 
tive. Though, as for the mean-field case, the magnitude 
of the polarization is conserved; fluctuations have renor- 
malized the quartic terms and shifted the phase bound- 
ary in Fig. [DJb). This behavior can also be recovered 
by a rcnormalization group analysis Q . One experimen- 
tal probe of the mctaclcctric transition is the changing 
behavior of the heat capacity C = -Td'^F/dT^. Before 
the metaelectric transition (small negative r) the relevant 
optic mode is "soft" and so the heat capacity follows the 
familiar Debye form C ~ T'^, whereas after the metaelec- 
tric transition (large negative r) , the relevant optic modes 
are "stiff" and so the heat capacity has an exponential 
dependence on temperature. At high temperature, in 
both cases the heat capacity has the expected classical 
behavior C = Sfce. 

Having conflrmcd the existence of a possible metaelec- 
tric behavior, wc now turn to consider the stability of the 
phase in the vicinity of the transition. Recent studies of 
itinerant ferromagnetism have suggested that such flrst 
order behavior can be preempted by the development of 
textured magnetic order analogous to that seen in the 
FFLO phase of superconductors This leaves open 

the question as to whether a textured phase can develop 
in the vicinity of the metaelectric transition. Our strat- 
egy to explore this possibility is to assume that the in- 
homogeneous phase is formed continuously, which allows 
us to develop a Landau expansion in the polarization <f> 
and texture wave vector Q. The onset of an inhomo- 
geneous phase is signaled by the coefficient of the ^^Q^ 
term turning negative. In our analysis we search primar- 
ily in the vicinity of the metaelectric transition at ^ = 
and consider a trial state with uniform polarization (f> 
that is for simplicity superimposed by an inhomogcneous 
component (f>cos(Q • r)(l, 1, 1). We then expand the free 
energy to quartic order in Q and discover that the pres- 
ence of a textured phase makes a contribution to the total 
energy of Q^^^l-u^^^ /6T:^^+7u^<i>^Q^ /60tt^^^]. Short 
of the first order transition where ^ < 0, the coefficient of 
is positive so the phase is not modulated. After the 
first order transition ^ turns positive driving the coeffi- 
cient of negative, revealing a finite Q instability in the 
region highlighted in Fig. [2Kb). The modulation carries 
polarization $ = r/2{u + 3/v). Though the analysis is 
restricted to the consideration of a potential continuous 
transition into the textured phase, and a simple form 
for the texture, it is sufficient to validate its existence. 
Refinements to include a putative first order transition 
or further textured phases would only enlarge the region 
of the phase diagram over which inhomogcneitics could 
be observed. Leaving aside potential textured phases we 



now turn to consider the behavior of the susceptibility 
across the phase diagram. 



B. Inverse susceptibility 

The inverse susceptibility provides an experimental 
window [l6l . [iTj onto the quantum critical properties 
of fcrroelectrics. Deep in the paraclcctric regime where 
i? = 7' + g/3 ^ (7d, the contribution to the inverse sus- 



ceptibility is X = d'^F/d4>^ 



R 



1(7- 



tan ^ 7) coth(^^), which is consistent with Barrett's for- 
mula [U for a gapped system. In the quantum critical 
regime we see three characteristic types of behavior for 
the inverse susceptibility 
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where the first term is from the diagonal contributions 
to Eqn. ([ll, and the latter two terms are the off-diagonal 
contribution, and ^ = qo/R. At low temperature ef- 
fects of long-range dipole interactions prevail as the off- 
diagonal fluctuating contribution rcnormalizcs the on- 
diagonal terms, with the linear temperature dependence 
of the term proportional to g giving a positive slope to the 
inverse susceptibility whereas the term could pro- 
vide a negative slope. At higher temperatures the con- 
tribution from the mean-field term dominates, which is 
also characteristic of quantum critical behavior and is in 
good agreement with recent experimental results Q . We 
note that the behavior is recovered by other models, 
including a diagrammatic resummation 0, 0] , the quan- 
tum spherical model [2^ . rcnormalization group stud- 
ies [13, [131, ^ self-consistent phonon model and an 
analogy to the temporal Casimir effect J9 L The behavior 
has also been observed experimentally |3l.ll6l. \v?\ . In both 
SrTiOs and KTaOa the initial linear negative slope and 
the quadratic ~ T"^ term conspire to cause a charac- 
teristic minimum in the inverse susceptibility. Using esti- 
mates for the parameters in Table [H the minimum occurs 
at T IK in both SrTiOa and KTaOa which is in good 
agreement with the experimental values of T = 1.6K and 
T = 3. OK respectively Finally, at high temperatures 
a classical term ~ T from the longitudinal fluctuat- 
ing term dominates from ^ lOOK which is again in good 
agreement with the experimental observations [^. 
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IV. DISCUSSION 

In this paper we have found that the polar fluctuat- 
ing phonons can drive a displacive ferroelectric through 
a first order metaelectric transition. Long range dipolar 
interactions did not alTect this critical phase behavior Q . 
However, long-range dipole interactions introduced into 
the action through the term (g — hq'^)(fP were pivotal 
in creating the correction to the inverse susceptibility 
~ —T that could explain the characteristic inverse 
susceptibility minimum |3|] , as well as provide important 
corrections to the self-consistent phonon treatment 0]. 

However, another mechanism, coupling of the soft op- 
tic modes to acoustic phonons could be significant. It has 
already been understood [1, 0, Q that a coupling with 
the acoustic phonons (p of the form —r]{y(p)(fP leads to 
a correction in of that could explain the char- 
acteristic minimum in the inverse susceptibility, and also 
has the capability of driving a first order transition [1, Q . 
This work and the results presented here motivate fur- 
ther experimental investigations into the inverse suscep- 
tibility and putative metaelectric transition that could 
shed light on the origin of the phase structure. Though 



the coupling to acoustic phonons complicated the solid 
state system, ultracold atoms in an optical lattice with 
long-range dipole interactions [l^l present a clean system 
that could provide a powerful tools to help unravel the 
properties of the generic Hamiltonian. 

One important simplification was to model the fer- 
roelectric with undamped dynamics. Damping would 
primarily arise due to free electrons, which can be in- 
troduced controUably through doping. Analogous to 
"avoided criticality" at a magnetic critical point which 
leads to non-Fermi liquid behavior and superconductiv- 
ity, ferroelectrics might also be expected to adopt novel 
behavior; for example doped SrTiOa [1^1, whereas un- 
dopcd SnTe [23| and GeTe [1^ become superconducting 
at low temperatures. This area presents a promising av- 
enue of research. Further open questions are to determine 
whether with just a change of parameters [1,[2^ the same 
formalism be applied to order-disorder ferroelectrics, and 
to consider the consequences of the coupling of fluctuat- 
ing polarization and magnetization that could arise in 
EuTiOa 

We thank Mark Dean, Gil Lonzarich, Stephen Rowley, 
and Montu Saxena for useful discussions. 
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